Transport through a quantum ring, a dot and a barrier 
embedded in a nanowire in magnetic field 
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We investigate the transport through a quantum ring, a dot and a barrier embedded in a nanowire 
in a homogeneous perpendicular magnetic field. To be able to treat scattering potentials of finite 
extent in magnetic field we use a mixed momentum-coordinate representation to obtain an integral 
equation for the multiband scattering matrix. For a large embedded quantum ring we are able 
to obtain Aharanov-Bohm type of oscillations with superimposed narrow resonances caused by 
interaction with quasi-bound states in the ring. We also employ scattering matrix approach to 
calculate the conductance through a semi-extended barrier or well in the wire. The numerical 
implementations we resort to in order to describe the cases of weak and intermediate magnetic field 
allow us to produce high resolution maps of the "near field" scattering wave functions, which are 
used to shed light on the underlying scattering processes. 



I. INTRODUCTION 

The influence of a single impurity on the conduc- 
tance of a quasi one-dimensional quantum channel has 
been investigated by several groups theoreticalljiSi 3 ^ 
and experimentally^ Commonly the impurities are con- 
sidered to be short range and represented by a ^-function, 
though treatments of more extended scatterers, like 
square barriers^ can be found. Recently, the applica- 
tion to nanosized systems, has spurred the use of general 
methods built on the Lippmann-Schwinger equation or 
the equivalent T-matrix formalism to describe the scat- 
tering of more general extended potentials in quantum 
channels 6 or curved wires. 7 

The inclusion of a constant homogeneous magnetic 
field perpendicular to the quasi one-dimensional electron 
channel or wire drastically changes the properties of the 
system. Without the magnetic field a centered symmetric 
scattering potential leads to "selection rules" that restrict 
the possible scattering processes. These are lifted by the 
magnetic field resulting in a rich structure contrasted to 
the conductance steps in an ideal wire as long as the mag- 
netic length is not much shorter than the width of the 
wire and the range of the scattering potential. 8 

The character of the Lorentz force does not allow us to 
establish a simple multimodc formulation of the scatter- 
ing process in configuration spaced but in strong mag- 
netic field Gurvitz used a scheme to develop a multi- 
mode formalism using a Fourier Transform with respect 
to the transport direction, and a truncation to a two- 
mode formalism allowed him to seek analytical solutions 
for a short range scatterer present in a wire with general 
confinement^ 

Here we extend this formalism by noting that in the 
case of a parabolic shape of the wire confinement we ob- 
tain coupled Lippmann-Schwinger equations with a non- 
local scattering potential in Fourier space for the different 



modes. A transformation of this system of equations to 
corresponding equations for the T-matrix shows that it 
bears a strong resemblance to the corresponding equa- 
tions for the system in no external magnetic field. 6 We 
exploit this fact to seek numerical solutions for the sys- 
tem in weak and intermediate strength of the magnetic 
field where a two mode approximation is not always war- 
ranted. One benefit of the numerical approach is that 
it allows us to map out with high resolution the proba- 
bility density for the scattering states near the scatterer. 
These "near field" solutions give us a good indication of 
the scattering process itself. We explore a quantum wire 
with an embedded quantum dot or a ring. We are able to 
increase the size of the ring to the limit that we observe 
Aharanov-Bohm type of oscillations. 

In order to investigate the scattering of potentials that 
are homogeneous in the direction perpendicular to the 
wire representing a barrier or a well we employ an al- 
ternative method based on mode matching. The smooth 
scattering potential is sliced into a series of (5-potentials. 
The scattering matrix is then constructed from repeated 
mode matching at each slice. 



II. MODELS 

We consider electron transport along a parabolically 
confined quantum wire parallel to the x-axis and perpen- 
dicular to a homogeneous magnetic field B = Bz. In the 
center of the wire the electrons are scattered by a poten- 
tial V sc (r) to be specified below. The system under inves- 
tigation is described by the Ffamiltonian H = Hn + V BC {^) 
with 
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where the wire is assumed to be parabolically confined, 
namely V c (y) — vri*^y 2 j2 with m* and f2 being, respec- 
tively, the effective mass of an electron in GaAs-based 
material and the confinement parameter. Here — e is the 
charge of an electron. We present two quite different ap- 
proaches to describe the scattering process of incoming 
electron states, one using wave function matching appro- 
priate to describe scattering from a thin homogeneous 
barrier or well perpendicular the wire, and the other one 
using a T-matrix formalism to calculate the scattering of 
an embedded quantum dot or ring in the wire. 



confinement parameter f2o* = U/(m*ao 2 ), and the mag- 
netic field B* — hc/(ea 2 ). As such, by defining the unit 
of the cyclotron frequency to be uj* = f^, the cyclotron 
frequency simply has the dimcnsionless form uj c = B. 

In the absence of scatterers the eigenfunctions of H 
can be written asm 



tp ± (x,y,k n ) 



x (y,k n ) 



(3) 



with ±fc„ being the wave vector along ±x in the nth 
transverse subband, where x^(y> k n ) satisfies the reduced 
dimcnsionless equation 



A. Scattering matrix approach 

In this section, we employ the scattering matrix ap- 
proach to the calculation of coherent electronic transport 
in a quantum wire in the presence of a Gaussian-profile 
scattering potential. The potential under investigation 
can be realized as a finger gate atop the wire and can be 
modeled by 



y sc (r) = l/ exp (-f3x 2 ) 
as is shown in Fig. ^ 
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X ± {y,k») = EnX^fakn), (4) 



which is a harmonic oscillator of frequency Q, w with a 
shifted center a n = w c fc n /fi^. These eigenmodes of the 
electron in a state described by ip(x, y, k n ) in a pure quan- 
tum wire have the energy spectrum 



E — E„ 



/C(fc„), 



(5) 



X (y, k n ) = N„H n V fl w (y =F a n ) 



composed of Landau levels E n = (n + 1/2) Cl w , where Q w 
= \/ uj 2 + SIq , shifted by the confining potential, and the 
kinetic energy /C(fc n ) = k 2 (rio/fl w ) 2 . The eigenfunctions 
X ± (y^ k n ) of the eigenmodes are given by 

exp — — (y T 0L n ) 
(6)' 

where N n — (Qw/ir) 1 ^ 4 (2 n n\) is a normalization 
constant. 

To utilize the scattering matrix approach, we divide 
the Gaussian scattering potential into a series of slices 
of width SL, each of them is described by a delta-profile 
potential 



V'sc (xi) = V SL exp {-(3x 2 ) 5[x- x{) 



(7) 



The Gaussian potential Eq. can thus be described 

by V sc (r) = YliHi ^sc {xi) if we divide the potential into 
sufficiently many pieces. For a right-going incident wave 
ip{x,y,k n ) from the nth. mode of the left reservoir, the 
corresponding scattering wave function can be expressed 
in the form 



ip$(x,y,k n ) = e lk " x xt(y,k n ) 
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FIG. 1: Barrier or well embedded in a quantum wire with Vo = 
(a) —6 meV and (b) 6 meV, respectively. Other parameters 
are B — T, a m = 33.7 nm, Eq = hQo = 6.0 meV, and 



1.897. 



To obtain a dimensionless expression we employ the 
Bohr radius do, and have thus the relevant units for the 



for x < Xi, and 



^\x,y,k n ) = V4 n e- ife «'* X + (y,fc„0 (9) 



for x > Xi. Following similar procedure we can also ob- 
tain the reflection and transmission coefficients, f n , n and 
t n i n , for a left-going incident wave. 
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The scattering should satisfy two boundary conditions, 
one requirement is that tp(x,y) has to be continuous at 
x — Xi and the other one stems from integration of the 
Schrodinger equation across x — Xi. Multiplying these 
boundary conditions by the eigenfunction x m {y) of the 
ordinary unshifted harmonic oscillator with confining fre- 
quency Cl w , and integrating over y, we obtain 

Imri (knO 4'n ~ ^mn' (&n') r n'n = ^mn (&n) (10) 
n' n' 

^ Jmn' i^n') Kl'n + ^mn' ^n') r n'n = ^"^rnn (^n) 
n' n' 

where the matrix elements are related to the overlap in- 
tegrals 

iL'= fxm(y)Xn>(y)dy, (12) 



expressed in terms of the incident electron energy E in 
the simple form 

O 2 N 

G(E) = —J2 T n(E), (17) 

n=0 

where N denotes the highest propagating mode incident 
from the source electrode. The current transmission coef- 
ficient T n (E) for an electron incident in the nth subband 
from the source electrode is given by 

t ^= E ^ra 2 - ( is ) 

n'(/i n />0) " 

The current reflection coefficient R n (E) can be calculated 
by a similar form to get current conservation condition 
for checking the numerical accuracy. 

B. T-matrix formalism 



J+ (fc n ,) = k n ,I+ n , (M + iV Q V+ n , (13) 



with 



and 



C = / Xm(y)v BC (y)Xn>(y)dy, (u) 



^mn' (&n') — kn<I mn i {k n ') ■ 



(15) 



Using Eqs. I|l() |l -I|ll |l and the corresponding equations 
for r l n , n and t l n , nl one can establish the scattering ma- 
trix for the Gaussian-shape potential.- From the ratio 
of the transmitted and the incident current we obtain 
the current transmission Tp a , in which a and /3 are, re- 
spectively, the incident and the transmitting lead. Since 
the current flowing into the terminal R (drain) from the 
system should take into account the current into the ter- 
minal L (source) from the system, the net charge current 
into terminal R is therefore given by 
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where f(E) indicates the Fermi distribution function and 
fi is the zero-bias Fermi level. The two-terminal con- 
ductance is given by G — 81/ dV. Below we assume 
that the scattering potential is located at the center of 
the nanowire, the applied perpendicular magnetic field 
is homogeneous, and the source-drain bias is sufficiently 
low. Hence, the zero temperature conductance can be 



We consider a quantum dot or a ring embedded in 
the wire and parameterize the scattering potential ac- 
cordingly combining two Gaussian functions of different 
shapes 



V sc (r) = Vi exp (-/3ir 2 ) + V 2 exp {-for 2 



(19) 



as is shown in Fig. El Together the magnetic field and the 
parabolic confinement define a natura l length scale a w — 
y/h/ (m*)tt w , where il w = y/w 2 + SIq, with the cyclotron 
frequency u> c = eB/(m*c), is the natural frequency of the 
quantum wire in a magnetic field. 

Along the lines of Gurvitz we choose to use the 
mixed momentum-coordinate presentation of the wave 
functions 9 



*b(p, y) 



dxip E (x,y)e lpx , 



and expand them in channel modes 

®e(p, y) = ^2 VnipJ&nip, y), 



(20) 



(21) 



i.e. in terms of the eigenfunctions for the pure paraboli- 
cally confined wire in magnetic field 



with the center coordinate y$ = pa 2 w (jj c /Vl w . These eigen- 
modes of the pure quantum wire have the energy spec- 
trum E np = E n + K n (p) with E n — KQ w (n + 1/2) and 



«n(p) 



( P a w ) 2 (nn ) 2 

2 



(23) 



Using Eq.'s (|20I23|I and performing a Fourier trans- 
form with respect to the coordinate x transforms the 
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where we have defined the effective band momentum 
k n (E) as 



(E-E n 



[k n (E)} 2 (HQq) 
2 ML,, 



(28) 



In the absence of a magnetic field it is possible to de- 
rive an equivalent effective one-dimensional multi-band 
Schrodinger equation to l|27l) in coordination spaced This 
multi-band equation is then usually transformed into a 
system of effective one-dimensional coupled Lippmann- 
Schwinger integral equations that is convenient for nu- 
merical computation. Here we can proceed along these 
lines, but the magnetic field forces us to do this in mo- 
mentum space where we shall see that the corresponding 
Lippmann-Schwinger equations are better transformed 
into integral equations for the T-matrix in order to fa- 
cilitate numerical evaluation. Considering Eq. I|27l) it is 
clear that the incoming scattering states satisfy 

[-(qa w ) 2 + (k n (E)a w ) 2 ] <£(q) = 0, (29) 

which implies a Greens function 

[-(qa w f + (k n (E)a w ) 2 } G n E {q) = 1. (30) 

The Greens function can now be used to write down cou- 
pled Lippmann-Schwinger equations in momentum space 



FIG. 2: Quantum ring (upper) or dot (lower) embedded in 
a quantum wire, B — T, a w = 33.7 nm, Eq — hQo = 1.0 
meV, /3%a% = 3.41, and /32a™ = 11-37. 



Schrodinger equation corresponding to the Hamiltonian 
l|T| into a system of coupled nonlocal integral equations 
in momentum space 

f dp 

Kn{q)<Pn(q)+^2 / 2^ V nn'(q,P)<Pn>(p) = {E - E n )lf n (q) , 
n' 

(24) 



where 



V nn '(q,p)= / dy(j)* n (q,y)V(q-p,y)<j> n ,{p,y), (25) 



and 



V(q -p,y)= dx e-^- p >V sc (x, v)- (26) 



The matrix elements (I25l2fi|) for the scattering potential 
(|19JI can be evaluated analytically since they consist of 
Gaussians and Hermite polynomials, see Appendix lAl 

The special form of the part of the energy dispersion 
K n (q) for parabolic confinement allows us now to rewrite 
Eq. <E3 as 



[-(qa w ) 2 + (k n (E)a w ) 2 ]ip n (q) 
2MI 



(Ml 



~]2 / ■^ V nn>(q,P)<Pn'(p), 

) ty-t I 



(27) 



cp n (q) = <p n (q)+G n E (q)J2 J 



dpa v 
~2V 



■V nn '(q,p)<Pn'(p), 



(31) 

where V nn >(q,p) = V nn '(q,p)2hQ w /[a w (tiQ. ) 2 ]. These 
equations are inconvenient for numerical evaluation as 
the in-state (p^ is proportional to a Dirac delta func- 
tion. Symbolically Eq. (|31|l can be expressed as ip = 
tp + GVip, and an iteration of the equation gives ip — 
<p° + GVtp°+GVGVip° + --- = (1 +GT)ip°, where we have 
introduced the T-matrix satisfying the symbolic equation 
T = V + VGT . Fully written the equation determining 
the T-matrix is 



T nn >(q,p) =V nn '(q,p) 



E 



dka u 
~2T 



■V nm >{q,k)G% (k)f nm ,(k,p). 



(32) 



This set of equations is easier to solve numerically than 
the equivalent Lippmann-Schwinger equations (|31(l after 
the singularities of the Greens function have been han- 
dled with special careii We obtain analytically the con- 
tribution of the poles of the Greens function and perform 
the remaining principal part integration by removing the 
singularity by a subtraction of a zeroii2ii^ 

Comparison with the nonseparable two-dimensional 
Lippmann-Schwinger equation in configuration space for 
extended scattering potential in a magnetic field gives 
the connection between the T-matrix and the probability 
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amplitude for transmission in mode n with momentum 
k„ if the in-state is in mode m with momentum k rn 



The conductance is then according to the Landauer- 
Biittiker formalism defined as 



G(E) = ^Tr[tt(£)t(£)], 



(34) 



where t is evaluated at the Fermi energy. 

Symbolically the wave function can be expressed as 
<p = (1 + GT)(p Q if the in-state tp° is given. Together 
with Eq. (J20J and J2U this gives 



iI>e{x,v) 



E 



<t>n{k n ,y) 

dqa 



2tt 



(35) 



for an incident electron with energy 25, in mode n with 
momentum fc„ . To calculate the wave function the same 
methods are used to isolate the contribution from the 
poles of the Greens function as were used for the calcu- 
lation of f with Eq. (|32J). 



III. RESULTS 
A. Embedded barrier 

In this section, we present our numerical results of ex- 
ploring electronic transport properties using a Gaussian- 
shape potential model described by Eq. (J2J - the con- 
ductance versus the incident electron energy E. The pa- 
rameters used to obtain our numerical results are taken 
from the GaAs-Al x Gai_ x As heterostructure system. The 
values that we choose for our material parameters are 
-ERyd = 5.93 meV and ao = 9.79 nm. 

The conductance Eq. I|17|l of the wire is presented 
in Fig. |3 for several values of magnetic field. To ex- 
plore the transport properties it is convenient to show 
the conductance as a function of energy of the incom- 
ing electron state scaled by the subband energy level 
spacing X = E/Ml w + 1/2 such that the integral val- 
ues of X indicates the number of incident modes. In 
Fig. |3 we present the conductance for magnetic fields 
with strengths from to 2.4 T for either weak (Vq = —6 
meV) or strong (Vq = — 12 meV) attractive potentials, 
as shown in Fig. [3Ja) and (b) , respectively. For the case 
of weak attractive potential shown in Fig.|3|[a), one can 
see that the dip structures in G(E) are pinned at around 
X ~ n + 0.85, and the location is insensitive to the mag- 
netic field. It turns out that these structures correspond 
to the electrons incident from subband n scattered elas- 
tically into the n + 1 subband threshold forming quasi- 
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FIG. 3: Conductance of a parabolically confined wire as a 
function of the energy parameter X = E/hQ w + 1/2 for 
various applied perpendicular magnetic fields. The ampli- 
tudes of the attractive impurity potential are (a) —6 meV 
and (b) —12 meV. Other parameters are hu>o = 6 meV and 
/3 2 a 2 w (B = 0) = 1.897. 



bound states^ It can be demonstrated that these quasi- 
bound states are formed in the leads out of the embedded 
Gaussian potential^ 

For the case of strong attractive potential shown in 
Fig. Otb) , one can see that there are two types of quasi- 
bound state features. The mechanism of sharp dips below 
the subband threshold is similar to the case of weak at- 
tractive potential. On the other hand, it is interesting to 
see the valley structures in Fig.|21for B ^ 0. These valleys 
correspond to quasi-bound states formed in the attrac- 
tive Gaussian potential. When the applied magnetic field 
becomes stronger, the blue shift of these valleys indicates 
such quasi-bound states are formed closer to the edge of 
the Gaussian potential due to the cyclotron motion. The 
large width of these valley structures implies the short 
life time of these quasi-bound states. When increasing 
the strength of the magnetic field, these valleys become 
wider. This indicates that the electrons in high mag- 
netic field with short cyclotron radius easily escape from 
the quasi-bound states formed in such a strong attrac- 
tive potential. We note in passing that in the absence of 
magnetic field, the intersubband transition is forbidden 
since the attractive potential is uniform in the transverse 
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direction, and we cannot see any dip structures in G(E). 
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FIG. 4: Conductance of a parabolically confined wire as a 
function of incident electron energy for various amplitudes of 
attractive potential. The other parameters are taken to be 
B = 1 T, Tiluo = 6 meV, and /3 2 a 2 w (B = 0) = 1.897. 

In Fig. 01 we study how the conductance can be af- 
fected by changing the amplitude of the attractive po- 
tential by fixing the strength of the magnetic field B = 1 
T, the confining potential huj y = 6 meV, and the Gaus- 
sian parameter /3 2 a^(-E> = 0) = 1.897. In the absence of 
the Gaussian potential (red curve), the electron transport 
manifests an ideal quantized conductance, the magnetic 
field plays no role. When increasing the amplitude of the 
attractive potential, the subband levels in the potential 
will decrease in energy. Therefore, we can find a red shift 
of the quasi-bound states. More precisely, for the cases 
of Vq = —3, —6, —12, and —18 meV, the dip structures 
occur at around E/hfl w = 1.95, 1.84, 1.54, and 1.17, 
respectively, in the attractive potential. It is also inter- 
esting to note that when the attractive potential is very 
strong, such as Vq — —18 me V, one can see a second dip 
structure appearing below the subband thresholds, both 
are quasi-bound states of the attractive potential. 

Figure |3] shows the conductance as a function of inci- 
dent electron energy for several values of magnetic field 
in the presence of a repulsive potential. The magnetic 
fields are tuned from to 2.4 T for either weak (Vq = 6 
meV) or strong (Vo = 12 meV) repulsive potentials, as 
shown in Fig. HJa) and (b), respectively. For the case 
of weak repulsive potential shown in Fig. EJa) , one can 
see that the conductance plateaus are suppressed from 
the ideal case. When increasing the applied perpendicu- 
lar magnetic field, the suppressed conductance plateaus 
tend to be enhanced back to the ideal case. For the case 
of strong repulsive potential (see Fig.[5fb)), the conduc- 
tance curves are suppressed much more then those for the 
weak repulsive one. Without a magnetic field, the con- 
ductance behaves according to Ohm's law. When the ap- 
plied magnetic field is increased, the conductance can be 
enhanced, and the quantization behavior becomes slowly 
recognizable. 

One may want to see how the conductance can be af- 
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FIG. 5: Conductance as a function of incident electron energy 
with various applied magnetic fields. The amplitudes of the 
repulsive potential barrier are Vo — (a) 6 meV and (b) 12 
meV. Other parameters are hu>o = 6 meV. and (3 2 a-w(B = 0) 
= 1.897. 
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FIG. 6: Conductance as a function of incident electron en- 
ergy with various amplitudes of repulsive potential. Other 
parameters are taken to be B = 1 T, Tiujo — 6 meV, and 
p 2 a 2 w {B = 0) = 1.897. 



fected by changing the amplitude of the repulsive poten- 
tial, as is shown in Fig. |5| We fix the strength of the 
magnetic field B = 1 T, the confining potential hujy = 6 
meV of the nanowire, and the potential configuration 
ho^uiB = 0) = 1.897. In the absence of scattering 
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potential, the conductance is ideally quantized. When 
increasing the amplitude of the repulsive potential, the 
electrons are reflected by the barrier and the conduc- 
tance is strongly suppressed. For weak barrier, such as 
Vq = 3 and 6 meV, the quantization behavior in G(E) is 
still visible. However, when the barrier is strong enough, 
such as Vq = 12, and 18 meV, not only the conductance 
is strongly suppressed, but the quantization behavior in 
G(E) has vanished. 

To conclude this section, we note in passing that when 
the scattering potential (well or barrier) is uniform in the 
transverse direction it does not break the translational in- 
variance along the lateral confining direction. However, 
in the presence of magnetic field, if such a scattering po- 
tential is a well then one can find quasi-bound states due 
to elastic intersubband transitions to a higher subband 
threshold. However, if the scattering potential is a bar- 
rier, one finds no quasi-bound state features even in a 
magnetic field up to 2.4 T. 



B. Embedded quantum ring and dot 

To model an embedded quantum ring with the param- 
eterization i|19|) we initially choose the parameters used 
in Fig. |2 such that when B — then /3x<z^ = 3.41, 
P2<it = 11-37, and %Q, W — 1.0 meV. (The parameters of 
the potential iJTSJ, (i\ and (3 2 do not depend on B, but 
a w does). V\ — —12 meV, and V% — 18 meV. We are 
thus investigating relatively broad wire with a small em- 
bedded ring structure with diameter of approximately 40 
nm. We assume the wire to be a GaAs wire as mentioned 
above. The conductance (|34|) of the wire is presented in 
Fig. |7| for several values of the magnetic field. To com- 
pare the results for various values of the magnetic field it 
is convenient to observe the conductance as function of 
the energy of the incoming electron state scaled by the 
distance of the energy subbands, i.e. E/(Ml w ) — E/E w 
and furthermore use X = E/Ml w + 1/2 such that an in- 
tegral value of X indicates the number of incident modes. 

In Fig. we see that as soon as the magnetic field is 
different from zero a strong Fano-likc 15-16 resonance dip 
appears in the first plateau just above X — 1.5. As we 
argue below the dip corresponds to a destructive quan- 
tum interference between a quasi-bound state in the ring 
and an in-state of the wire. Figure |S] displays the total 
probability to find an electron in the wire close to the 
scattering center, the quantum ring. The probability is 
calculated using the wave function (|35|l for two values of 
the energy of the incoming electron in the lowest trans- 
verse mode, n = 0. Just below the resonance at X = 1.4 
Fig. |BK reveals us a normal scattering process. The scat- 
tering only takes place very close to x ~ and on the left 
hand side we see the interference pattern for the incoming 
and the reflected wave. On the right hand side the elec- 
trons only travel in one transverse mode and only to the 
right so we have a constant probability already a short 




B=1.6 T 
B=2.2 T 
B=2.8T 



2 3 
X 



FIG. 7: The conductance of a parabolic wire with an em- 
bedded ring in units of Go = 2e 2 /h. E w — h£l w , Vi = —12 
meV, V 2 = 18 meV, Mlo = 1.0 meV, f3ia 2 w (B = 0) = 3.412, 
Paat(B = 0) = 11.37, and 9 subbands are included. 



distance away from the scattering center. The situation 
is quite different in Fig. \j^p that displays the probability 
density for the state exactly in the resonance dip. Here 
no transmitted wave is present, but the probability close 
to the quantum ring is high enough that the probability 
for the incoming and the reflected waves is not visible on 
the color scale used. 

The symmetry of the quasi-bound state indicates that 
it is an evanescent state belonging to the second subband 
n = 1. Without a magnetic field the scattering via the 
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FIG. 8: The probability density of the scattering state 
ipE{x,y) in the parabolic quantum wire in the presence of 
an embedded quantum ring (Fig. corresponding to the 
conductance in Fig. |7Ji at B = 0.1 T. The incident energy 
X — 1.4 (a), and X = 1.538 (b) corresponding to the dip in 
the conductance. 



evanescent state in the second subband is forbidden in 
the case of a symmetric potential placed in the middle of 
the wire. In that case a dip occurs in the second band 
due to a scattering through a evanescent state in the third 
subbandi 2 i 6 i 17 

In order to further support our view that the resonance 
is due to a quasi-bound state of the quantum ring located 
in the continuum of the first subband (The ring lowers 
a state in the second subband into the first one), we see 
in Fig. how the broadening or narrowing of the wire 
has little effects on the energy of the state. On the other 
hand, as seen in Fig. ^jp the energy of the quasi-bound 
state changes linearly with the depth of the ring poten- 
tial. 

For some intermediate values of the magnetic field we 
see other minima occurring in the conductance closer to 
the end of the first step. For example, for B — 0.8 T this 
is visible in Fig. |3> at X = 1.933 and 1.991. The cor- 
responding probability densities are seen in Fig. 1101 The 
symmetry of both densities indicates that the dips are 
caused by scattering via evanescent states of the second 
subband just as the dip in the middle of the first conduc- 
tance step. These states are quasi-bound states of the 
ring further in the continuum of the first subband. The 
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FIG. 9: The conductance of a parabolic wire with an embed- 
ded ring in units of Go = 2e 2 /h as a function of the wire 
confinement Eo — hQo (a), and depth of the ring Vi (b). 
V 2 = 18 meV, pia? w (B = 0) = 3.412, /3 2 a 2 w (B = 0) = 11.37, 
and 9 subbands are included. E w — KQ, W , B — 0.8 T and 
Vi = -12 meV in (a). E = 1.0 meV and B = 0.1 T in (b). 

higher state, Fig.llOb has acquired more of the character 
of the geometry of the wire than the ring, and it extends 
far beyond the ring. The presence of Fano lineshapes 
in the conductance is not surprising as the mesoscopic 
Fano effect was already experimentally reported for both 
a single electron transistor— and an Aharanov-Bohm in- 
terferometer with an embedded quantum dotii2*22i Nev- 
ertheless, in these two experiments the wire was much 
smaller than the mesoscopic system that caused the Fano 
interference. The results presented here suggest that this 
may be observed also in the case of a broad wire. 

At still higher magnetic field, see Fig. 0;, the conduc- 
tance has approached the ideal case as the magnetic field 
has now squeezed the wave functions together and closer 
to the edge as soon as the momentum is different from 
zero. The wave function thus bypasses the scattering po- 
tential. We shall see this effect clearer below. 

The conductance of a wire with an embedded dot is 
presented in Fig. ^] The effects of an increasing mag- 
netic field become very clear if we compare the probabil- 
ity density for the dips at X = 1.679 when B = 0.6 T, 
and the one at X = 1.557 when B = 1.2 T, see Fig. Ill 
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FIG. 10: The probability density of the scattering state 
ipE{x,y) in the parabolic quantum wire in the presence of 
an embedded quantum ring (Fig. corresponding to the 
conductance in Fig. 0> at B = 0.8 T. The incident energy 
X — 1.933 (a), and X — 1.991 (b), corresponding to two 
minima in the conductance at the end of the first step. 



B=1.2 T 
B=1.4T 
B=1.6T 



(b) 



FIG. 11: The conductance of a parabolic wire with an em- 
bedded dot in units of Go = 2e 2 /h. E w — TiQ w , Vi = 12 
meV, V 2 = -18 meV, titto = 1.0 meV, f3ial(B = 0) = 3.412, 
f3 2 al{B = 0) = 11.37, and 9 subbands are included. 



Both cases show a partial blocking of the channel due to 
backscattering caused by a quasi-bound state created by 
an evanescent state of the second subband, but the main 
difference is the total separation of the incoming and the 
reflected channel at the higher magnetic field. At the 
lower magnetic field we still see a interference pattern 
between these channels. 

To explore further the scattering from the embedded 
dot when there is not a complete separation between the 
edge and bulk channels we show the probability density 
for X = 2.235 at B = 0.6 in Fig. [B] This energy cor- 
responds to a dip seen in Fig. Illb . Due to the magnetic 
field and the scattering potential there is always a scat- 
tering between these two channels irrespective of whether 
the instate belongs to the n — (Fig. 113b ) or the n = 1 
(Fig. 113b - ) mode. This is visible in the probability density 
with interference pattern in all channels. 

The situation is completely different at the higher mag- 
netic field B = 1.2 T seen in Fig. Here the same 
scaled energy as before, X — 2.235, corresponds to a 
peak in the conductance displayed in Fig. lllb . The value 
of the conductance peak indicates that there is very lit- 
tle backscattering. The edge channel (n = 0) is almost 
entirely tranmitted as Fig. 114b shows, but a quasi-bound 



state is seen in Fig. 114b belonging to the same subband 
as the instate. 

The small quantum ring embedded in the broad wire 
(Fig. |2J) is to small to show any indication of Aharanov- 
Bohm oscillations, and as the magnetic length gets 
smaller with increasing field strength the edge states by- 
pass the ring. To change this situation we also did cal- 
culations for a larger ring shown in Fig. 1151 compared to 
the smaller ring used in the preceeding calculations. 

Of course the parabolic confinement of the wire always 
leads to the situation that at a high enough magnetic 
field the edge states will not be scattered by the quantum 
ring potential, but now at an intermediate field strength 
the magnetic length compares more favorably with the 
size scale of the ring as can be seen in the conductance 
displayed in Fig. [TB] for both B = and 1 T. At B = 1 
T we see oscillations growing in wavelength with E or 
k n (E) both for mode n — and 1. 

The oscillations in the conductance at B = 1 T are 
caused by a simple geometrical resonance where the 
wavelength of the scattering state in the ring has to com- 
pare appropriately with the circumference of the ring 
to build constructive or destructive interference, i. e. a 
Aharanov-Bohm like effect. This also explains the grow- 
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FIG. 12: The probability density of the scattering state 
ipE{x,y) in the parabolic quantum wire in the presence of an 
embedded quantum dot (Fig.[3Jl, corresponding to the conduc- 
tance in Fig. ^3 The incident energy X = 1.679 at B — 0.6 
T (a), and X = 1.557 at B = 1.2 T (b), corresponding to two 
minima in the conductance at the end of the first step. 



FIG. 13: The probability density of the scattering state 
tl>E{x,y) in the parabolic quantum wire in the presence of an 
embedded quantum dot (Fig. corresponding to the con- 
ductance in Fig. Illfa at B — 0.6 T for a state with incident 
energy X — 2.235, in mode n — (a), and n = 1 (b). 



ing wavelength of the oscillation with (E — Eg). Even 
though the same condition lies at the root of the energy 
spectrum of stationary states in a ring in equilibrium we 
are not probing here the energy spectrum of the ring. We 
would like to mention that a similar oscillatory behavior 
of the conductance as a function of the energy was re- 
ported by Sivan et al. 21 in the case of a quantum dot in 
high magnetic field. Since in high magnetic field the cy- 
clotron radius is smaller than the ring radius one expects 
the electrons to travel within the ring along skipping or- 
bits before they leave through the wire. 

The large size of the embedded ring in this case and 
it finite depth mean that quasi-bound states will not be 
of the same simple structure as seen for the smaller ring. 
This can be verified by the probability densities shown 
in Fig. E|f° r the two dips at X = 1.319 and 1.347, and 
for the maximum at 1.425. The total transmission of the 
only mode, n = 0, in Fig. 117b causes the perfect left right 
symmetry, but the probability density in Fig. 117b corre- 
sponding to the dip at X = 1.319 reflects the asymmetry 
caused by the confining parabolic potential to the ring 
seen in Fig. 115b . The structure of the evanescent states 
in Fig. 117b and b indicates that they are caused by states 
in the third and fifth energy bands and probably also 



states in higher bands. This persistence of eigenstates 
or scarring of wave functions in open systems has been 
discussed by Akis et &t~ for quantum dots, and here we 
confirm it for an open quantum ring. 

Superimposed on the Aharanov-Bohm like oscillations 
in the conductance in Fig. 1161 we have narrow resonances 
that are caused by interaction with quasi-bound states 
of the ring. In Fig. 118b we show the probability density 
for the scattering state for a Aharanov-Bohm peak at 
X = 1.46. The density for a minimum in the oscillation 
is similar except for the addition of the reflected wave. In 
Fig. 118b we display the probability density for the first 
narrow resonance seen, at X — 1.135. Here we can iden- 
tify a long lived evanescent state in the second subband. 



IV. SUMMARY 

We have successfully extended a multiband transport 
formalism build on the Lippmann-Schwinger equation in 
a magnetic field to be able to describe an unbiased trans- 
port through a broad wire with embedded small or large 
quantum dots and rings defined by a smooth potential. 
The calculation of the probability density for the scatter- 
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FIG. 14: The probability density of the scattering state 
ipE{x,y) in the parabolic quantum wire in the presence of an 
embedded quantum dot (Fig. corresponding to the con- 
ductance in Fig. Illb at B = 1.2 T for a state with incident 
energy X — 2.235, in mode n — (a), and n = 1 (b). 
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FIG. 15: A contour plot of the potential of an embedded 
quantum ring in wire, (a) a ring with f3ia%, = 3.41, and 
P 2 at = 11.37, corresponding to the ring in the upper panel 
of Fig. (b) a large ring with /3ial = 0.0682(B=0), and 
P 2 al(B = 0) = 0.682. E = hQ, w = 1.0 meV, a w = 33.7 nm 
at B = T. 



ing states allows us to shed light on internal processes and 
resonances that in some cases reflect interaction between 
states in several subbands of the wire. We observe well 
known evanescent states and Fano resonances produced 
by these interactions. In the case of a large ring with 
finite width we observe Ahranov-Bohm type of oscilla- 
tions superimposed with narrow resonances reflecting its 
energy spectrum. 

Due to the wide range of system parameters used we 
had to pay extra attention to the accuracy of the numer- 
ical methods employed. 



APPENDIX A: MATRIX ELEMENTS OF THE 
SCATTERING POTENTIAL 

The matrix elements of the potential represented by a 
single Gaussian function V = Vq exp (— (3 x x 2 — f3 y y 2 ) is 



according to Eq.'s (|25I26|I 



V nn '(q,p) = V a w \/ — e:q> 



{q-pf 



Inn'(q,P), (Al) 



where 



I nn ,(q,p)= dyct>* n (q,y)e-^ 2 ct> n ,(p,y). (A2) 




FIG. 16: The conductance of a parabolic wire with a large 
embedded ring (corresponding to Fig. I15b ~l in units of Go = 



2e 2 /h. E w = hQ w , Vi = -12 meV, V 2 
meV, /3ial(B = 0) = 0.0682, p 2 ai(B 
subbands are included. 



is meV, nn = 1-0 

0) = 0.682, and 13 
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FIG. 18: The probability density of the scattering state 
ipE{x,y) in the parabolic quantum wire in the presence of an 
embedded large quantum ring (Fig. 115b ) . corresponding to 
the conductance in Fig. 1161 at B — 1 T for a state with inci- 
dent energy X — 1.46 corresponding to a peak (a), X = 1.135 
in a narrow dip (b). The incoming mode is n = 0. 



Insertion of the expressions for the eigenfunctions IL'l'l) 
yields 



FIG. 17: The probability density of the scattering state 
ipE{x,y) in the parabolic quantum wire in the presence of an 
embedded large quantum ring (Fig. 115b ). corresponding to 
the conductance in Fig. 1161 at B — T for a state with inci- 
dent energy X = 1.319 corresponding to a dip (a), X = 1.347 
in a dip (b), and X = 1.425 at a maximum (c). The incoming 
mode is n = 0. 



Inn'{q,p) 



exp 



■(s n + s n ,) 2 _ s 2 n +s 2 n , 



4C 



2 n + n W Cn\n'\ 
EE (p) ^H p (-V2s n )H p (-V2s n ,) 

min(n— p,n —q) 

E *« 

(=0 




(A3) 



where N = n + n' - p - q, s„ = yfi = k n a w uj c /fl w , 
z = (s n + s n >)/{2VC), C = (1 + p y al), and b = (1 - 
2/C). When the variable & assumes negative values the 
combination (v / fc) JV_2 '-ffjv-2/(' ■ ■ /Vb) still supplies the 
correct real value. 
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